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ABSTRACT 

The baryon acoustic oscillation (BAO) feature in the two-point correlation function of galaxies 
supplies a standard ruler to probe the expansion history of the Universe. The current surveys 
covering a comoving volume sufficient to unveil the BAO scale are limited to a redshift z < 1 . 
We study here several galaxy selection schemes, aiming at building an emission-line-galaxy 
(ELG) sample in the redshift range 0.6 < z < 1.7, that would be suitable for future BAO 
studies, providing a highly-biased galaxy sample. We analyze the angular galaxy clustering 
of galaxy selections at the redshifts 0.5, 0.7, 0.8, 1, and 1.2 and we combine this analysis 
with a halo occupation distribution (HOD) model to derive the properties of the haloes these 
galaxies inhabit, in particular the galaxy bias on large scales. We also perform a weak lensing 
analysis (aperture statistics) to extract the galaxy bias and the cross-correlation coefficient and 
compare the HOD prediction. 

We apply this analysis on a data set composed of the photometry of the deep co-addition 
on SDSS Stripe 82 (225 deg^), of CFHT/Stripe82 deep / band weak lensing survey and of the 
WISE infrared photometric band Wl. 

The analysis on the SDSS-III/CMASS selection at z = 0.5 is in agreement with previous 
studies on the tracer, moreover we measure its cross-correlation coefficient r - 1. 16+0.35. For 
the higher redshift bins, we confirm the trends that the brightest galaxy populations selected 
are strongly biased {b > 1.5), but we are limited by current data-sets depth to derive precise 
values of the galaxy bias. A survey using such tracers of the mass field will guarantee a high 
significance detection of the BAO. 

Key words: cosmology - large scale structure - galaxy - bias - baryonic acoustic oscillations 
- dark matter 
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1 INTRODUCTION 

Working within the standard cosmological framework described by 
the Friedmann model, recent analyses of cosmic microwave back- 
ground (CMB) and type la supemovae (Snia) data confirm the pic- 
ture of our Universe as almost spatially flat, with an energy content 
known at the percent level and dominated by dark components, and 
undergoing a phase of accelerated expansion since the last 7 Gyr 
{e.g. [Komatsu et al.|2011[ [Suzuki et al.|2012) . A confirmation of 
these results is provided by the use of the baryonic acoustic oscil- 
lations (BAO) signature in the clustering of galaxies ( ,Frieman et al.[ 
[2008[ >, which offers a standard ruler for the measurement of the 
angular distance-redshift relation. Current measurements of BAO 
are limited to redshift z < 0.7 and at precision larger than 2%. 
[Beutler et al.[ ( |20l"T| l provides the latest measurement at z = 0.106 
on the 6dF Galaxy Survey ( [Jones et al.|2009 i with a precision of 
4.5%. [Blake et al.[j201 1] > and [Anderson et al.[ ( |20T2} made the lat- 
est measurements at z 0.6 using, the WiggleZ Dark Energy Sur- 
vey (Drink water et al.|2010 ) and the Sloan Digital Sky Survey III 
Baryonic Oscillation Spectroscopic Survey (SDSS-III/BOSS here- 
after) jEisenstein et al.[[2011[ [Dawson et al.[[20T3| l, respectively. 
The detection at the sub-percent level of BAO at larger redshift, 
0-7 ^ z ^ 2, constitutes the next step toward the investigation of 
the accelerated expansion of the Universe. 

The measurement of BAO in the power spectrum of galaxies 
relies on the assumption that 'galaxies are tracers of the underly- 
ing matter distribution' l |Kaiser|1984^ . The validity and the limits 
of this assumption are not yet completely understood theoretically 
or observationally. To constrain cosmology to the sub percent level 
with the measurement of the BAO feature in the galaxy clustering, 
it is necessary to understand the relation between the galaxy distri- 
bution and the dark matter field. 

Although the BAO scale is in the linear regime of structure for- 
mation ( [Peebles & Yu[1970[[Eisenstein & Hu[1998[ >, the systematic 
errors on its measurement are dominated by non-linear structure 
growth at small scales and by redshift-space distortions. These fea- 
tures broaden and shift the peak in the two-point correlation func- 
tion, which is now a standard method to measure galaxy clustering 
( [Eisenstein et al.[2007^ . Once the properties of the BAO tracers are 
known, it is possible to reduce these systematic errors below 0.2%; 
see [Padmanabhan & Whit e (2009 ); Mehta et al. (201 1 1. Therefore, 
it is of prime interest to investigate, on small and medium scales, 
the clustering properties of future BAO tracers and their relation 
to dark matter. Currently, to bridge the gap between dark matter 
haloes and galaxies, there are different methods. 

The first approach is to use a model that provides how the 
galaxies are distributed in haloes (Halo Occupation Distribution 
model; hereafter HO D); see [Dekel & Lahav] ([T999ll; [Ma & Fry[ 
( [2000) ; [Seljak[ ( |2000| l; [Co"oray & Sheth[ j2Q02 ). The HOD model 
is fitted on the measure of the galaxy clustering to predict the prop- 
erties of the haloes that the galaxies inhabit. In particular, HOD pre- 
dicts the galaxy bias; see [White et al.[p011[ i; [Zehavi et al.1 ( [2011[ i; 
[Parkinson et al.[ ( [2012a^ for the HOD applied to galaxies used in 
BAO surveys. 

The alternative method is based on weak lensing (WL). With 
recent surveys, it has become possible to combine on small areas (a 
few deg^) different statistical measures of the galaxy distribution: 
the galaxy clustering, the cosmic shear, and the galaxy-galaxy lens- 



are being developed to interpret the combination of both measures 
( [Leauth aud et a l.|2011[|Cacciato et al.|2012[[More et al.|2012]|van[ 
[den Bosch et al.|2012| >. 

The application of such analysis on thousands of square de- 
grees by combining weak lensing with galaxy clustering is however 
currently limited by the depth and the point spread function (PSF) 
characterization of the photometry ( [Huff et al.|2011[ l. 

Another possibility to determine the galaxy bias is the direct 
comparison with large simulations of dark matter using halo abun- 
dance matching technique (HAM), ( Kravtsov et al.|2004"{ [Vale &[ 



[Ostriker|2004[[Conroy et al.|2006[[Trujillo-Gomez et al.|201 1\ . Ap- 
plied to the analysis of WiggleZ and SDSS-III/BOSS observations, 
this technique shows similar results as the HOD; see [Prada et al.[ 
( |20T2t ; |Nuza et al.[ ( |20T2l l; [Parkinson etaL] ( [2012b| >. 

In this paper, we investigate the clustering amplitude of galax- 
ies selected to be the BAO tracers of future large spectroscopic 
surveys as SDSS-IV/eBOSS, BigBOSS, DESspec, PFS-SuMiRe, 
Euclid ( [Comparat et a l. 2013, Schlegel et al. 20 lT[|Laureijs et aTj 



201 1 1. The clustering amplitude being proportional to oc (fccrg) , we 



consider a constant erg = 0.81 (value from CMB analysis WMAP7 
[Komatsu et al.[ ( [201 1[ >) in order to derive the value of the bias b. 
For this analysis, we combine an HOD approach with a weak lens- 
ing analysis. We obtain rather quantitative rough estimates of the 
galaxy bias, that should not be taken at face value. We make the 
point for a high clustering amplitude (or a high bias), justifying tar- 
get selection scheme of future BAO surveys and a more in depth 
investigation to make forecasts for these surveys accounting for bi- 
ases at these levels. In fact signal-to-noise for BAO studies scales as 
the density of tracers times the clustering amplitude of these trac- 
ers, and hence BAO studies profit from highly clustered tracers. 

In section|2]we examine the galaxy bias and the two different 
methods used to measure it: the halo occupation distribution model 
and the weak lensing aperture statistics method. We also discuss 
the limits of each measurement. In section [3] we describe the data 
on which our measurement is based. In section|4]we discuss the bias 
measurement from clustering and weak lensing for each tracer and 
the associated errors. Finally, in section[5]we discuss the feasibility 
of future BAO studies using the described tracers. 

Throughout this paper we assume a flat ACDM cosmology 
characterized by WMAP7 parameters Hq = 100 km s"' Mpc"', 
with h=0.71, OcDM = 0.226, flj, = 0.0455, n^ = 0.7285, erg = 
0.81, "initial = 0. 966 (|Komatsu et al.[20Il^ . Magnitudes are given in 
the AB system ( [Oke & Gunn|1983| l. 



ing. Using COSMOS data, (Capak e t al.|2007||Scovil le et al. 2007 
[Ilbert et al.[2009| l [Jullo etal.| ( |2012^ ; |Leauthaud et al.| ( |2012 ) led the 
way in combining these measures to constrain the galaxy-halo re- 
lation and the evolution of the galaxy bias. Improved HOD models 



2 GALAXY CLUSTERING - WEAK-LENSING JOINT 
ANALYSIS 

Galaxy clustering and gravitational lensing probe, respectively, the 
galaxy overdensity field Sg(r,z), which is a discrete random vari- 
able function of position r and redshift z, and the matter overden- 
sity field, (5m (r, z), which is a discrete random variable related to the 
latter by some functional form generically called "bias", eventually 
allowing fo r stochasticity ([T egmark & Peebles|I998[ [Dekel & La-[ 
hav 1999, Tegmark & Bromley 1999)[jAs random variables, only 
the A'-point correlation functions (or functions thereof) are mean- 
ingful and can be measured; by combining galaxy clustering and 
gravitational lensing one can measure the auto-correlation func- 
tions {6f) and (Sl^) and the cross-correlation (dgd„,) and eventually 



' For clarity, we drop in the following the space and time dependence from 
6s. 
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define the (linear) bias, bg, and stociiasticity, r, parameters, after 
IHoekstra et ari ( |2502l l 



(SI) ' 



{6g6m) 



(1) 



A positive (negative) value for r indicates correlation (anticorrela- 
tion) between the galaxy and the matter fields, while r = l(r = -l) 
indicates total correlation (total anticorrelation). 



2.1 bg by clustering, an HOD approach 

To derive the galaxy bias, we need the auto-correlation of matter 
and of galaxies. The galaxy auto-correlation function (5^) is esti- 
mated with |Landy & Szalay| ( |1993) 'minimum variance estimator'. 
Then we use a HOD to interpret the measure of auto-correlation, as 
it provides an insight into the relation between the galaxies and the 
haloes they inhabit. 

The main ingredients of the model are n(M, z) the halo mass 
function i.e. the density of haloes of mass M present at redshift 
z\ the bias of haloes compared to matter bi,iM,z)\ and the auto- 
correlation of the matter. The parameter of the model is a func- 
tion that gives the number of galaxies present in a halo of mass 
M, designated N(M). For a given N(M) the HOD model outputs 
the corresponding auto-correlation function for galaxies. A fit on 
the galaxy auto-correlation function is performed to derive the an- 
alytical form of N{M). During the fit, the cosmological parameters 
are fixed at their fiducial value (WMAP7). The large scale galaxy 
bias of a volume-limited sample of galaxies described by N{M) is 
obtained with equation (|2j. 

/ n(M,z)bh(M,z)N(M)dM 



j n(M, z)N(M)dM 



(2) 



For more details about the implementation of the HOD we use in 
this study, see |Cooray & Sheth| ( |2002) ; |Tinker et"^ ( (2005l l; [Coupon| 
|et al.| 2012^ . We compute correlation functions with the software 
and the HOD fits are performed using CosMoPMCQ see 



[Kilbinger et dl ( [MT| [20091 |20Tol i; [Coupon et al.| ( |20T2l i 

The main limit of the model is the input auto-correlation of 
matter. This information is taken from [Smith et al.| ([2003} and 
has errors of ~ 8% at small scales and < 3% at large scales 
k < lOh Mpc"' in the ACDM paradigm. Also the HOD model is not 
valid in different cosmological paradigms and can therefore not be 
used to rule out cosmologies: it can only be used to replace galaxy 
samples in the galaxy clustering history, which is the aim of this 
paper. 



galaxies contains the information of the cross-correlation (SgS,,,). 
The auto-correlation (6^) of galaxies is derived in the same way as 
for the HOD. 

To be able to compare the auto- and cross-correlation re- 
sults, the measurements are convolved with Bessel kernels to trans- 
form the correlations in aperture correlations. This formalism is 
described in Schneide r et al.| ( fl998[ l and applied on RCS data in 
Hoekstra et al. ( 2002 ), GaBoDS data in Simon etaljl|2007}, and on 
COSMOS in Julio et al. (2012) . We use the same routines as in | Julio | 
jet al.| ( [2012^ . The first output is the galaxy bias as the ratio of the 
aperture auto-correlation of galaxies, denoted N, and that of mat- 
ter, denoted Map, in an aperture of an angular radius d; see equation 
l[3j. The second measurement is the cross-correlation coefficient r, 
that measures the randomness; see equation l[4l. 



fc(0) = /,((?, n,„,nA) 



(3) 



(4) 



The functions /[ and /2 are normalizations. They correct for the 
fact that different cosmological volumes are probed by the differ- 
ent statistics. The functions are computed assuming a jSmith et al.j 
( |2003^ non-linear power spectrum with unbiased foreground galax- 
ies (b - r = 1), that uses the |Eisenstein & Hu[ ( |1998| l transfer func- 
tion, constrained by the WMAP7 cosmological parameters ( |Ko-| 
[matsu et al.|2011t . 

The WL also uses the [Smith et~aL] 12003) 1 non-linear power 
spectrum and has therefore the same precision limitations as the 
HOD. Moreover, an additional constraint is that, to be able to ex- 
tract the WL signal, the sources must be behind the lenses. Because 
of this point, we are not able to have a significant measure of the 
matter galaxy cross-correlation by WL for the faintest galaxy sam- 
ples (see section[3|. 



3 DATA 

In this section, we present the data sets that are used for this anal- 
ysis. First, we explain the catalogs used to target galaxies tracing 
BAO. Secondly, we describe the catalogs containing the galaxy 
shape measurements. We then describe the dark matter halo simu- 
lation we used to check the consistency of the mass measurements 
from the HOD. Finally, we detail the selection functions for the 
BAO tracers. 



2.2 bg and r by weak-lensing 

Using weak-lensing analysis, we derive both bg and r. For this pur- 
pose, we need the three auto and cross-correlations mentioned be- 
fore. The WL measurement is based on the analysis of the distri- 
bution of shapes of the galaxies located behind the BAO tracers. 
Their intrinsic shapes are distorted by the foreground mass around 
the BAO tracers. The auto-correlation of the ellipticities of the 
background galaxies contains the information of the matter clus- 
tering between the background galaxies and us, {6j„). The cross- 
correlation of the ellipticities with the position of the foreground 



3.1 SDSS-Stripe 82 galaxy catalog 



http://www2.iap.fr/users/kilbinge/athena/ 
http://www2.iap.fr/users/kilbinge/CosmoPMC/ 



The Sloan Digital Sky Survey ( [York et al.|2000[ [Eisenstein et al.[ 
[201 1[ >, delivered under the data release 9 (DR9, |Ahn et al.|20T2l i7 
covers 14,555 square degrees in the 5 photometric bands u, g, r, 
i, z. It is the largest volume multi-color extragalactic photometric 
survey available today. The 3cr magnitude depths are: u = 22.0, 
g = 22.2, r = 22.2, / = 21.3; see[Fukugita et"aT] ( [T996} for the 
description of the fi lters, |Gunn et al. jl998 1 for the characteristics 



of the camera, and [Gunn et al.[ ( |2006 !) for the description of the 
2.5m telescope located at Apache Point Observatory (New Mexico, 
USA). 

The data set used here for the selection of the different BAO 
tracers is the Sloan Digital Sky Survey deep Stripe 82 coadd ( |Annis[ 
et al.|20lT I. It covers an area of 225 deg^ in one continuous stripe 
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Table 1. Samples of BAO tracers used in this analysis, 'b' stands for bright 
and 'f ' for faint, z is the mean redshift of the sample, <t- is the dispersion 
of the redshift distribution, N the total amount of tracers in the sample, Ns 
the density of tracers per deg^, '± scale' is the transverse scale in /r"' Mpc 
con'esponding to one degree at |. 



sample 


z 


O". 


N 


Ns 
deg"2 


± scale 
/i"' Mpc deg 


CMASS 


0.53 


0.10 


22k 


98 


16.0 


LRG-WISE b 


0.58 


0.12 


15k 


70 


17.8 


LRG-WISE f 


0.69 


0.13 


59k 


275 


18.5 


ELG gri b 


0.80 


0.08 


119k 


530 


19.2 


ELG gri f 


0.81 


0.15 


226k 


1007 


19.3 


ELG ugri b 


0.95 


0.24 


156k 


693 


20.3 


ELG ugri f 


0.94 


0.22 


436k 


1938 


20.2 


ELG ugr b 


1.28 


0.36 


77k 


341 


21.5 


ELG ugr f 


1.26 


0.34 


237k 


1052 


21.5 



10= 



10" 



10= 



backgrount^'^ourc^s ' 




CMASS 
LRG-WISE 

ELG gri 

-,LG ugri 

ELG ug 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1 



Figure 1. Redshift distributions of the samples. For each sample, solid 
(dashed) lines refer to the bright (faint) sample. 22 000 CMASS (red), 
15 000 (59 000) LRG-WISE (blue), 119 000 (226 000) gri ELG (green), 
156 000 (436 000) ugri (yellow), 77 000 (237 000) ugr ELG (magenta). 
The dashed black line represents the redshift distribution of the background 
sources used in the weak-lensing analysis. 



-1.25° < ^j2ooo < 1.25°, -34° < 0-12000 < 46° with five photometric 
bands u, g, r, i, z. We also use the infrared counterpart of this area 
observed by WISE at 3.4yum Wright et al.|2010[ >. Table[T] gives an 
overview of the tracers selected. 

The photometric redshift distribution of the BAO tracers are 
obtained using the CFHT-LS Wide W4 field photometric redshift 
catalog T0007. The photometric redshift accuracy is estimated to 
be cr- < 0.08(1 -Hz) for i 24. The data and cataloging methods are 
described in |Ilbert et al.]p006) ; |Coupon et al.| ( |2009^ the T0007 re- 
lease documenj^ The CFHT-LS W4 field overlaps the stripe 82 on 
a few square degrees. Both datasets have similar depths in magni- 
tude. From a position match of the targets selected on the stripe 82, 
we obtain a sub-catalog of the targets with photometric redshifts, 
from which we infer the redshift distribution of the selections. The 
redshift distributions are shown on Fig.[T] 

With current data, it is not feasible to obtain a reliable pho- 
tometric redshift on the complete stripe 82. Therefore, we cannot 
construct clean volume limited samples, as required by the HOD 
analysis. We can only make magnitude limited samples. The im- 
pact of the galaxies located in the tails of redshift distribution on 
the angular clustering is discussed in paragraph |4.1| 



http://www.cfht.hawaii.edu/Science/CFHTLS/ 



3.2 CFHT-Stripe 82 weak-lensing catalog 

The existing comprehensive data sets available on the SDSS Stripe- 
82 area were complemented with high-quality i'-band observations 
at CFHT. The main purpose of this one-band Megaprime@CFHT 
(see [Boulade et al.||2003 l follow-up survey is to enrich already 
existing multi-colour data with deep optical observations suit- 
able for weak lensing studies. The CFHT Stripe-82 program 
(CFHT/Stripe82) is a large collaborative effort between the Cana- 
dian/French and Brazilian communities!^ observations encom- 
pass 169 MegaPrime new pointings resulting in an effective survey 
area of about 150 deg^. This area comprises 169-4(due to fring- 
ing)=165 pointings CFHT/Stripe82 and 8 CFHT-LS Wide point- 
ings i.e. a total of 173 pointings. All the data were obtained under 
excellent seeing conditions, and the final co-added images show an 
image quality between 0.4 and 0.8 arcsec with a median of 0.59 
arcsec. Each pointing was obtained in four dithered observations 
with an exposure time of 410s, each resulting in a 5-cr limiting 
magnitude in a 2 arcsec diameter aperture of about 24. 

The data processing, which starts with the Elixir preprocessed 
images obtained from CADC and ends with co-added images and 
all necessary quantities to perform weak gravitational lensing stud- 
ies, follows closely the descriptions given in |Erben et al.| ( (20T2l l. 

We exploit the opportunity that the CFHTLS includes the HST 
COSMOS survey field, and verify the calibration of our shear mea- 
surement pipeline for ground-based data against an independent 
analysis of much higher resolution space-based data. We stack the 
subsets of the CFHTLS DEEP2 imaging to the same depth as the 
CFHT Stripe 82 survey and analyze it using KSB method ( |Kaiser| 
|et al.|1 995 1. The KSB method is used to measure the shapes of all 
the selected galaxi es. Our implementation of KSB is based on the 
KSBf9(|3 pipeline jHeymans et"aL]|2006b 1. This has been generi- 
cally tested on simulated images containing a known shear signal 
as part of the Shear TEsting Programme (STEP: [Heymans et al.| 
(2006a); M assey et al.| ( |2007^ ) and the GRavitational lEnsing Accu- 
racy Testing challenge (GREAT08: [Bridle et al.|pOTO) ; GREATIO: 
[Kitching et al.| l [2012^ ). The systematic errors of this pipeline are 
well controlled in all cases. This pipeline has also been used on the 
shear measurement of CFHTLS-Wide Wl fields ( |Shan et aL|20T2l ). 
Finally, we match galaxies in the CFHT/Stripe82 catalog to those 
of the HST COSMOS catalog ( |Leauthaud et al.|2010^ to calibrate 
the measurements. 

Here is the description of the calibration tests we made on 
the lensing catalog. The multiplicative calibration component is 
(mi, nil) = (0.13,0.14) and the additive component is (C[,C2) = 
(0.0003,0.0015). These coefficients are obtained from simulations 
and connect through the linear relation between the observed 
galaxy ellipticities and the true ones (gobs = (1 + m)[y + etrue] + c). 
They mainly depend on the algorithm that measures the shapes. In 
order to check for significant additive systematics we performed 
four tests. 

(i) We test the average shear (y) across all 173 pointings of the 
CFHT/Stripe82 field. Fig. |2] demonstrates that the average shear is 
almost consistent with zero as expected, for galaxies of all sizes, 
magnitudes and signal-to-noise ratios (SNR). As in |Heymans et al.| 
( |2012^ , (§2) is clearly dependent on the SNR (red points). We 
quantify the additive component of the calibration correction with 



5 The survey CFHT identification numbers are 10BB09, 10BC22, and 
10BF23 



http : //www . roe . ac . uk/ -heymans/KSBf 98/Home . html 
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Figure 2. Average shear measurements (gi) (left) and (§2) (right) of KSB 
from galaxies in ('-band observations of the entire CFHT/Stripe82 field. In 
the absence of additive systematics, these should be consistent with zero. In 
practice, they always remain within the dashed lines that indicate an order of 
magnitude lower than the 1-10% shear signal around clusters. Top, middle 
and bottom panels show trends as a function of galaxy magnitude, size and 
signal-to-noise ratio. The black points denote the average shear from all the 
galaxies in the two shear catalogs. The red points denote (§2) before the 
correction described in (i). 



a simple relation (^2) = a + bSNR. We find a best fit reduced 
= 0.54 with a = -0.002, = 0.000118. The corrected average 
shear is almost consistent with zero as expected, for galaxies of all 
sizes, magnitudes and signal-to-noise ratios (SNR). Furthermore, 
the average shear of each pointing has also been tested, see Fig. 
[3] We find that the mean shear (gi) of pointing S82p24p is around 
-0.0126 ± 0.0017, which is far from zero compared to other obser- 
vation pointings. Systematic error residual signal may exist in this 
pointing. We remove the related galaxies in the final catalogue. 

(ii) Star-galaxy cross correlation (residual systematics due to 
imperfect PSF). For the total 173 pointings, we find our PSF cor- 
rection is well within requirements for our analysis. Fig. |4] shows 
the correlation ^syj between the corrected shapes of galaxies and 
the uncorrected shapes of the stars. We normalize the star-galaxy 
ellipticity correlation by the uncorrected star-galaxy ellipticity cor- 
relation to assess its impact on shear measurements l [Bacon et al.| 
|2003[ l; see Fig. [5] Though, we find ^,,,5 is obviously higher than ^+ 
in one pointing S82p35p at scale 9 > 3arcmin, which reveals the 
presence of a systematic error residual signal in this pointing. We 
remove the related galaxies in the final catalogue. 

(iii) Two-point shear correlation functions. We calculate the 
two-point shear correlation functions of the shear catalogue and 
compare them to an analytical prediction. The comparison between 
the analytical prediction of halo model and the measured ^+ is good, 
see Fig.|6] Note the cross-correlation <f,x between the tangential and 
rotated shear should be zero. 

(iv) £,i?-mode signals. We test for the (nonphysical) 6-mode 
signal. The S-mode signal coiTesponds to the imaginary component 
of P^(£). We rotate all galaxy shears through 45° and remeasure the 



Figure 3. Average shear (gi) of 173 pointings, in which S82p24p field is 
biased from zero. 



Mr 
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Figure 4. The cross-correlation between shear measurements and stellar 
ellipticities as a function of the separation between galaxies and stars, aver- 
aged throughout the CFHTLS Stripe 82 field. The black and red points are 
shear correlation function and star-galaxy cross correlation, respectively. If 
all residual influence of the observational PSF has been successfully re- 
moved from the galaxy shape measurements, the red points should be con- 
sistent with zero. 



ii-mode signal. We measure that the i?-mode signals are around 
zero except for small scales; see Fig.|7] 

For the lensing catalog, we only select galaxies fainter than the 
foreground samples by 1 magnitude in i-band, so that we have 
no overlap between the two samples. The complete redshift dis- 
tribution of the sources used for the weak lensing analysis is dis- 
played on Fig.[T|with a thick dashed black line labeled 'background 
sources'. 
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Figure 5. Difference 6( between shear auto-correlation ^-i- and star-galaxy 
cross congelation ^sys at scale = 3arcmin. The difference delta of the 
S82p35p pointing is higher than others. 



10 

9[arcmin] 

Figure 7. The aperture mass from the entire CFHT/Stripe82 field. Red 
points show the B-mode, black points show the E-mode. The error bars 
of the E-mode include statistical noise added in quadrature to the non- 
Gaussian cosmic variance. Only statistical uncertainty contributes to the 
error budget for the B-mode. The errorbars are Icr confidence level. Within 
2(T errors the B-modes are consistent with zero. 




S[arcmin] 



Figure 6. The two-point shear correlation functions of the shear catalogue: 
f+ (black points), f_ (red points), and ftx (blue points). Solid lines show 
the analytical predictions of f + . Note the cross-correlation f tx between the 
tangential and rotated shear should be zero. 



3.3 Multidark-Stripe 82 halo catalog 

We use the Multidark dark matter simulation to characterize quali- 
tatively our samples to the first order using halo abundance match- 
ing (HAM). To perform the match, we use snapshots of the Multi- 
dark simulation corresponding to the mean redshift of our tracers. 
Every snapshot occupies a volume of 1 h^^Gpc^ filled with haloes 
obtained with bound density maximum algorithm with a cut-off at 



200 times the critical density (logm M2oo/Mo)- These data sets are 
taken from the database Multidarl^( Riebe et al.|201 1| >. 

In each snapshot, we select the dark matter halo population 
that corresponds to the number density of galaxies with z- cr- < 
z < z+cT; selected on the Stripe 82. This links each galaxy sample to 
a halo sample. The density, velocity cuts and characteristic masses 
and satellite fractions are given for each sample in Table |2] 

The HAM brings two components to this analysis. First it pro- 
vides an upper limit for the mass of the haloes containing the BAO 
tracers, which constitutes a good consistency check for the HOD 
models. Secondly, it provides the scatter of the mass distribution. 
Given that we measure clustering over ~ 200 deg^ and weak lens- 
ing on ~ 150 deg^ and not the full sky, we use this scatter to display 
an error on the mean halo mass deduced from the maximum likeli- 
hood HOD model. Thirdly, we use the satellite fraction of the halo 
samples as a lower limit to the satellite fraction that we should find 
for the associated galaxy sample. 

From the qualitative numbers obtained by the HAM, the LRG- 
WISE samples appear to be located in massive halos of ~ 10'^^ Mq 
with a low satellite fraction, similarly to the CMASS tracers. Con- 
cerning the ELG, the satellite fraction are higher and the halo 
masses upper limits are lower: < IO'^ ^Mq for the bright samples 
and < 10'""* Mq for the faint samples. 

3.4 BAO tracer selections 

The latest BAO tracers used at z > 0.5 are SDSS-III Constant Mass 
galaxies (CMASS) at z ^ 0.53 ^Dawson et al.|2013^ and W iggleZ 
Emission Line Galaxies at z - 0.6 Prinkwater et al.|20101 ). The 
next stage BAO experiments plan to use Luminous Red Galaxies 
(LRG) and Emission Line Galaxies (ELG) at z > 0.7 ( Schlegel] 
|et al.|20l"T|>. The target se lection of ELG tracers is extensively dis- 
cussed in jComparat et al.| ( |2013^ . 

^ http://www.multidark.org/ 
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Table 2. Abundance matching of the samples of BAO tracers used in this 
analysis, 'b' stands for bright and 'f ' for faint. The density of tracers is in 
units of 10""* ft^Mpc"^. Vd, or Vmax is the velocity cut applied to the mul- 
tidark simulation snapshot to match the density d. The obtained sample is 
characterized by the three quartiles of its distributions given in the columns 
logjQ M2oo/Wo- fsai is the fraction of satellites in the halo catalog. 



sample 


density 
d 


Vcir km s ' 


logioM2oo/Mo 


fsat 


CMASS 
LRG-WISE b 
LRG-WISE f 


2.64 
1.57 
5.06 


384 
431 
324 


13.17t34 

13.28!°;;^ 
1 7 Q9+0.36 


5% 
4% 

7% 


ELG gri b 
ELG gri f 
ELG ugri b 
ELG ugri f 
ELG ugr b 


8.16 
13.37 

5.44 
16.45 

1.45 


292 
251 
317 
238 
435 


,2 77+0.37 

12.59+«;^^ 
12.84!«;35 

13.21!«;?^ 


8% 
10% 
7% 
11% 

3% 


ELG ugr f 


5.76 


315 


12.81!«;J4 


7% 



The selection of the tracers used in this analysis is only based 
on photometric criteria, depth and color. We do not use the pho- 
tometric redshifts for the selection as their quality varies strongly 
across the stripe. The selection criteria target the most luminous 
galaxies in incremental redshift intervals to z = 2; see Fig.[8]which 
compares all samples in absolute magnitude M, versus redshift or 
rest-frame M„-Mg. The cleanest way to select BAO tracers would 
be to select by stellar mass, but with current data this is not possi- 
ble. 



3.4.1 CMASS at z =^ 0.5 

The CMASS tracer is the SDSS-m/BOSS experiment BAO tracer 
( [Dawson et al.|2013)l. Th e current BAO detection using the data re- 
lease 9 I Ahn et al.|2012^ (a third of the observation plan) with the 
CMASS tracers has a 6.7cr significance ( Anderson et al.|20T2l l. Its 
selection is described in detail in|Dawson et al. (2013). The parent 
catalog of CMASS selection on stripe 82 contains 22034 tracers 
i.e. ~ 98 deg"^. We use the complete CMASS selection, not only 
the galaxies confirmed by spectroscopy, in order to avoid fiber col- 
lision issues. The redshift distribution is shown by the red solid line 
in Fig[T] The mean redshift is 0.53 with a dispersion of 0.1. We use 
snapshot number 62 of the multidark simulation. The density of the 
tracer at z is 2.64 x 10"''/i'Mpc"^, which corresponds to a circular 
velocity cut of Vcii = 384.1 km s"'. This value of Vdr is consistent 
with the study done by Nuza et al.| (|2012'). The halo mass distribu- 



tion peaks at log|Q M2oo/A^o 
fraction expected is 5%. 



13 17+0-34 



The minimum satellite 



3.4.2 LRG-WISE at z ^ 0.65 

This selection is the extension in redshift of the red CMASS selec- 
tion. We name this sample the 'Luminous Red Galaxies selected 
with WISE' (LRG-WISE); see Schlegel et_ai](goTTji, Prakash et al. 
(2013), in preparation. On stripe 82, it selects 15643 galaxies i.e., a 
density of 70 deg"- at a mean redshift of 0.58 for the bright sample 
with 1 < 20 and ~ 59000 galaxies at a mean redshift of 0.69 for the 
faint sample with i < 21 i.e., a density of 262 deg"-. The selection 
is as follows : 



• SDSS z-band flag must be true, and WISE channel 1 flag must 
beO 

• Objects with extreme astrometric uncertainties were rejected 

• bright selection : (; < 20 or z < 20) and r — i > 1 and 
r-W\>2{r-i)- 0.5 and r > 16 and r < 22.5 

• faint selection : (i < 21 or z < 21) and r - i > 1 and r-W\> 
2(r - i) - 0.5 and r > 16 and r < 22.5 

The redshift distribution of the bright (faint) samples are shown 
with a solid (dashed) blue line on Fig.[T| Their densities are 1.57 
and 5.06 x 10^^/i'Mpc"', respectively. We use the multidark snap- 
shot 58 for the bright and the faint samples HAM. The samples 
correspond to cuts at V^i = 431.5 km s"' and 323.6 km s"', re- 
spectively. The halo mass distributions peak at log^ M2OO/M0 = 



13.28^]^, 



andlog,oM2oo/Mo 
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respectively. They con- 



stitute an upper limit for the halo mass these tracer inhabit. The 
minimum satellite fractions expected are 4% and 7%, respectively. 



3.4.3 £LG gri af z ^ 0.8 



The gri ELG selection jComparat et al.|[2013) is an intermediate 
selection between the reddest galaxies selected in SDSS-III and 
the bluest ones selected by the ugr and ugri schemes (see below), 
identifies emission-line galaxies with emission lines a small Balmer 
break. On Stripe 82, it selects a 119 000 and 226000 galaxies for 
the bright and faint samples, which correspond to a projected den- 
sity of about 530 and 1000 deg"- at mean redshift z - 0.8 and 
0.81, respectively. The bright sample is limited by 19 < ; < 21.5, 
and the faint sample by 19 < / < 22. The densities are respectively 
8.16 and 13.37 xlO"'* h^Mpc^^. We use multidark snapshot number 
56 to perform the HAM. The samples match with a circular veloc- 
ity cut of 292.2, 250.7 km s"'. The halo mass distributions peak 



at logio M200/M0 
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and logio M200/M0 = 12.34!; 



+0.44 

0.25' 



respectively. They constitute an upper limit for the halo mass these 
tracer inhabit. The minimum satellite fractions expected are 8% and 
10%, respectively. 



3.4.4 ELG ugr at z=^ 1.2 



The ugr ELG selection ( [Comparat et al.|2013| l selects the blue tail 
of the color distribution aiming for galaxies with strong emission 
lines at higher redshift (z - 1.2). The ugr redshift distribution is 
flatter and extends to higher redshifts. The bright sample is limited 
by 20 < g < 22.5, and the faint sample by 20 < g < 23. For 
the bright and faint sample it selects, 77 000 and 237 000 galax- 
ies, corresponding to a density of 341 and 1052 deg"^, respectively. 
Their mean redshift are at 1.28 and 1.26. We use multidark snap- 
shot number 48. The bright sample density, 1.45 x 10"''/i'Mpc"', 
corresponds to a cut at 434.5 km s"', which gives an upper limit 
of logioM2oo/Mo < 13.2l^°-]l. The faint sample density, 5.76 x 
10"''/j^Mpc"^, corresponds to a cut at 315.0 km s"', which gives an 
upper limit of log^ M2oo/Afo S; 12.81!55^g. The minimum satellite 
fractions expected are 3% and 7%, respectively. 



3.4.5 ELG ugri at z ^ 0.95: building SDSS-IV/eBOSS 

The ugri ELG selection is designed to be the BAO tracer sample 
for SDSS-IV/eBOSS project (eBOSS collaboration 2013 in prepa- 
ration). It is a mix of the two previous gri and ugr ELG selections. 
The selection criteria are 



we use dereddened r,i,z and 3.4pm (WISE 1) magnitudes 



20 < g < 22.5 for the bright; 20 < g < 23 for the faint 
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redshift redshift Mu-Mg Mu-Mg 

Figure 8. Absolute magnitude M, vs redshift and rest-frame M„-Mg. All galaxies with i < 23.5 (grey), CMASS (pink), LRG-WISE (light blue), gri ELG (red), 

ugri ELG (green), ugr ELG (blue). First panel shows the bright selections are quite luminous with M, 23. The second panel shows the faint selections 

are less luminous than the bright selections with M, 22. The third panel shows the bright selections span the complete rest-frame color range from ugr 

ELG (blue) to the CMASS (pink). The fourth panel shows the faint selections. Their rest frame colors are similar to that of the bright selections, except for the 
LRG-WISE faint galaxies that are bluer than the bright sample. The ugri ELG overlap with the gri and ugr ELG selections, which is expected as the selection 
encompasses galaxies from both other selections. This figure was made with the same CFHT W4 matched sample as for the redshift distributions. 



• -0.5 <ii-r< 0.7(g - -t- 0.1 and g-i> -0.5 

The mean redshifts are 0.95 and 0.94. The bright sample con- 
tains 156000 galaxies and the faint sample 436000 (i.e. densi- 
ties of 693 and 1938 deg -). We use multidark snapshot num- 
ber 52 for the HAM. The densities are respectively 5.44 and 
16.45 xlO"''/i''Mpc"', corresponding to mass upper limits of 
logioM2oo/Mo < 12.84!«;35 and log,o Mjoo/Me < 12.48!°;^3. The 
minimum satellite fractions expected are 7% and 11%, respectively. 



4 MEASUREMENTS 

In this section, we present the measurements made through cluster- 
ing and weak lensing analysis of the data. The correlation function 
measurements are displayed on Fig s. |9|13| Figs . 1 1 4| and 1 1 5 1 show 
the galaxy bias measurement. Table |4. 1 1 summarizes the details of 
the measurements. 



4.1 Angular clustering measure and HOD fits 

We measure the angular clustering, denoted w(d), for angles in the 
range 0.001° to 1.5°. In this manner we avoid the signal due to 
blended sources. In the case of LRG-WISE, the clustering analysis 
starts at 0.01°. We use between 15 and 20 angular bins regularly 
log-spaced. The errors and the covariance matrix are issued from 
112 bootstrap realizations of the angular clustering measure. Fig. 
|9] shows a comparison of the angular clustering of all previously 
defined samples. 

The most precise estimations of the angular correlation func- 
tions are between 0.01° and 0.1°, due to a balance between the den- 
sity of BAO tracers and the geometry of the survey. For the sparse 
tracers (the bright selections at low redshift), below 0.01°, there are 
only a few pairs of tracers. For the densest tracers, the correlation 
function is precise down to 0.001°. For all tracers above 0.1°, given 
the geometry of the survey (a long stripe), the correlation function 
measurements are covariant. We correct from the integral constraint 
as in Wa ll & Jenkins |pOT2t . 

The main trends in the angular clustering measurements fol- 
low expectations. The lower redshift samples are more clustered 
than higher redshift samples. The bright samples are more clustered 
than the faint samples in particular at small scales. 



To interpret the angular clustering, we use the HOD 
parametrization given in pi, where erf(x) = f ^ e^' dt is the 'er- 
ror function'. From N(M) we deduce the angle averaged 3D corre- 
lation function, that we project using the Limber equation jLimber| 



1954||Simon|2007 i to obtain whod- The^^ accounting for the ac- 



curacy of the fit is defined by equation ^ 
^ log M - log M, 



N(M) : 



1+erf 



1 + 



M-Mo\ 
Ml j 



(5) 



- (l^obscrvcd - >Vh0d)''^.C ' .(Wobserved - WhOd) + 



(A'; 



.gal 



(6) 



'Ngal 



The main source of error is the uncertainty on the number 
of galaxies representing the redshift distribution of each popula- 
tion. This uncertainty lies between 5 and 20% for the bright tracers 
(CMASS and ELG bright). For the faint tracers it is between 30 
and 40% (error denoted cr^gai). The second source of errors is the 
statistical errors on the measurement of w(d), which is accounted 
for by the covariance matrix C. Another source of errors is the cos- 
mic variance. The fractional error on the measured density due to 
cosmic variance is cr„ < 0.1; see [Newman & Davis| ( |20()2) . We 
account for this when computing the error on the HOD derived pa- 
rameters. We add to the HOD error the dispersion on the halo mass 



obtained with the HAM crj^j^f^y 



(Ml,alo)+CrHOn(A^halo) = (Mhalo). 



In practice, crHAM(^iiaio) » cThod^^^mo) we therefore only use 
'^HAM(*^haio); sec Table 



4.1 



We perform additional tests to search for systematic errors. We 
test if the angular clustering measured is representative of the red- 
shift distribution. We compute the angular clustering on only the 
tracers that have a reliable photoZ. We find a discrepancy on the 
angular clustering smaller by a factor of three than the error on the 
angular clustering measurement. We therefore assume the angular 
clustering measurement is representative of the redshift distribu- 
tion. 

The validity of the use of an HOD model can be discussed. 
N{M) is valid to predict the clustering when applied to a vol- 
ume limited sample. Current data does not allow us to construct 
a clean volume limited sample, as we do not have good photomet- 
ric redshifts for all the tracers. We cannot cut a clean sample with 
z 6 z ± CTj. We quantify the impact of the wide redshift distribu- 
tion on the angular clustering on large scale, where it matter for the 
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determination of the galaxy bias, using the Multidark simulation. 
We compare the angular clustering of haloes HAM selected with a 
top-hat redshift distribution (1 if z 6 z±it-, else 0) and the observed 
redshift distribution. The large scale angular clustering 6 > 0.08° is 
unchanged. The small scale angular clustering does change. There- 
fore we assume the large scale angular clustering measured on the 
full sample is representative of the sample located in z ± cr., and 
that the bias, within the error bars will be correct. The values of 
the galaxy bias derived by the HOD method are given in Table |4. l] 
and plotted in Figs. |14|and|15[ This shows the major flaw of this 
HOD parametrization, that cannot account for a galaxy population 
with a wide redshift distribution. A hint to try and model this would 
be to introduce a redshift dependence for the parameters describing 
N{M). We further discuss this issue in section[5] 

4.2 Weak lensing measure 

We measure N{ff)Map(d), N^(9), M^pid) with the same routines as 
in |Jullo et al.| ( |2012) . The signal-to-noise ratio obtained on the bias 
is ~ 3 for tracers with z < 1 and of ~ 1.5 for z > 1. Due to this 
low signal to noise measurement, we exclude the ugr ELG sample 
at z ^ 1.2 from the WL analysis. The measurements are shown on 
Figs.[T0](CMASS),[TT](LRG-WISE),[T2](^n- ELG),[T3](Mgn ELG). 
We derive the galaxy bias for all these samples; see Table |4.l] The 
cross-correlation N(d)M„p{6) is measured with precision only for 
the lower redshift CMASS sample. Therefore r is computed only 
for this sample. 

The main source of error arises from the error on the shape 
measurement; see the data description section. The uncertainty on 
the redshift distribution can also potentially introduce some system- 
atic errors. All statistical errors are propagated with bootstrap as de- 
scribed in |Jullo et al.| ( |2012) l. For each foreground and background 
galaxy catalog, 100 samples are randomly drawn and processed to 
compute the aperture auto and cross-correlation functions as well 
as the bias and correlation coefficient. The background source cata- 
logs are cut at 1 magnitude fainter than the foreground. Systematic 
error on Ml^iG) on large scale due to the lack of information at 
small scale is below 1% ( [Kilbinger et al.|2006^ . 

4.3 Cross-correlation coefficient and galaxy bias 

Measurements of the galaxy clustering as a function of the lumi- 
nosity were performed in the local universe by ( [Coil et al.|2008| 



'Skibba & Sheth"2009' "Zehavi et al."2011;'Guo et al."2012V They 
showed the clustering amplitude is greater for more luminous sam- 
ples. Thus we expect to measure high values of the galaxy bias for 
all samples. There is no straightforward existing relation that can 
predict the expected galaxy bias. 

4.3.1 Cross-correlation coefficient 

We measure the cross-correlation coefficient only for the CMASS 
sample, r = 1.16 ± 0.35. This measurement is a first, and implies 
the CMASS galaxies are fully correlated with the matter field. This 
strongly supports current BAO analysis made on this sample \An-\ 
[derson et al.|2012| l. 

With current data, it is not possible to constrain the cross- 
correlation coefficient of the other galaxy samples. 

4.3.2 Galaxy bias 

For the CMASS, LRG-WISE bright, gri ELG bright, and ugri ELG 
bright samples, the measurements of the two methods are in agree- 
ment; see Fig. [14] We find the bias is above 1.5 for the all ELG 
selections, which is consistent with the measurements made by 
[Coupon et al.| ( |2012| l on its brightest samples. We consider both 
measurements give a trend about the size of dark matter haloes 
these galaxies inhabit; see Fig. [15] Comparing the values of the 
bias with the curves of the bias for a constant halo mass after ITin-1 
|ker et al.l ( |2005| l and with the mass distributions obtained via HAM 
shows all these samples are tracing the most massive haloes at their 
respective redshifts. 

For the LRG-WISE faint, gri ELG faint and ugri ELG faint 
samples, there is a systematic disagreement between the methods : 
the clustering analysis outputs a higher bias than the weak lensing 
analysis; see Fig.[T4j The values are consistent at Icr. Thus for the 
LRG-WISE faint, gri ELG faint and ugri ELG faint samples, we 
find the combination of current surveys "CFHT/Stripe82 + CFHT- 
LS + SDSS-S82 deep coadd" is too shallow to derive a reliable 
value of the bias. We need a fair spectroscopic sample of such trac- 
ers to improve the photometric redshift distribution and reduce the 
error on the redshift distribution and a deeper survey for the shape 
of background galaxies. 

For the ugr ELG bright (and faint), with current data, the WL 
analysis is not feasible. The galaxy bias derived using only the clus- 
tering information is although of prime interest. It is consistent with 
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Table 3. Results from the clustering and the weak lensing analysis. The number of degrees of freedom is ~ 18 for the HOD fits. The values of the^^ show 
current HOD models fail in explaining angular clustering of the ELGs. The column < Mj states the relation M;,o/„ — Mj. fodusterine is the galaxy bias obtained 
with the HOD, (Z)wl) with the weak-lensing, and (^wl) is cross-correlation coefficient. 



sample maximum likelihood HOD parameter ^^/dof halo mass ^'clustering (^wl) (''wl) 
logi'lS^) log(^) ^Qg(B^) (^logM a (Mi,ah) < Mj 



CMASS 


13.433 


14.576 


10.73 


0.242 


1.790 


0.3012 


13 


779 


> 


2.5 ± 0.5 


2.05 ± 0.3 


LRG-WISE bright 






















1.77 ±0.28 


LRG-WISE faint 




















2.0 ±0.2 


1.44 ±0.22 


ELG gri bright 


13.51 


14.51 


11.07 


0.89 


0.81 


16.05 


13 


219 




1.95±0.4 


2.06 ± 0.39 


ELG gri faint 


13.256 


14.931 


11.994 


0.268 


0.691 


3.379 


13 


558 


> 




1.57 ±0.38 


ELG ugri bright 


12.41 


12.043 


12.422 


0.58 


0.64 


2.88 


13 


157 




2 0+0.31 


1.73 ±0.33 


ELG ugri faint 


13.049 


14.443 


13.089 


0.69 


0.54 


6.95 


13 


084 


> 




1.56 ±0.34 


ELG ugr bright 


13.179 


13.409 


13.074 


0.63 


0.62 


1.54 


13 


241 


< 


9 7+0.3 
^■'-0.2 




ELG ugr faint 


13.103 


14.386 


12.928 


0.592 


0.779 


6.53 


13 


126 









1.16 ±0.35 




9[deg] »[deg] 9 [deg] 

Figure 10. WL measurements for the BOSS-CMASS sample. From the left to the right the galaxy auto-correlation N^(d), the matter autocorrelation Ml^iO), 
and the galaxy-matter cross-coirelation N(0)Map(8). The black line represent the mean value and the red lines the Icr error contours. Measurements of N'(6) 
and (d) are clean between 0.05° and 1°. N(9)Map(d) is clean between 0.02° and 0.2°. We thus can deduce bg and r. 
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9[deg] »[deg] 9 [deg] 

Figure 12. WL measurements for the gri ELG bright (top row) faint (bottoin row) sample. From the left to the right the galaxy auto-correlation N^(ff), the 
inatter autocorrelation Mlp(6), and the galaxy-matter cross-correlation N(6)Map(6). The black line represent the inean value and the red lines the Icr eiTor 
contours. Measurements are clean between 0.1° and 1° for the M^p(6) and N^(8) of the bright sample and N'(6) of the faint sainple. The other lensing 
measurements are not robust. 




9[deg] 9[deg] 9 [deg] 

Figure 13. WL measurements for the ugri ELG bright (top row) faint (bottom row) sample. From the left to the right the galaxy auto-coiTelation N^{6), the 
matter autocorrelation M^ip(6), and the galaxy-matter cross-correlation N(6)Map(6). The black line represent the mean value and the red lines the Icr error 
contours. Measureinents are clean between 0.05° and 1° for the M^p(6) and N'(6) of the bright sample and N^(6) of the faint sainple. The other lensing 
measurements are not robust. 



the fact that the ugr ELG bright sample occupies the most massive it is always 0.5 dex higher than the upper limit set by the HAM. 
halos at z ~ 1.2. The clustering curves of the brightest galaxies in each redshift bin 

therefore constitute a new challenge to halo occupation models. 



4.3.3 Satellite fraction, halo mass and HOD 

The shape of the angular clustering for the ELG samples suggests a 
high satellite fraction. Indeed there is a strong bump for 9 < 0.02°. 
The HOD models fail to reproduce trustworthy satellite fractions. A 
similar effect is seen for the mean halo mass predicted by the HOD; 



5 DISCUSSION 

The measurements presented in this paper are meaningful only 
for the CMASS, LRG-WISE bright, gri ELG bright, ugri ELG 
bright, and ugr ELG bright (clustering only) samples. The follow- 
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Figure 15. Large scale galaxy bias of the different BAO tracers vs. redshift. The interquartile are shown as error bars (25% - 75% intervals). The dashed black 
lines correspond to the bias of haloes of mass log[Q M200/M0 of 12.5, 13 and 13.5 from bottom to the top, after |Tinker et al.|j2005^ . The figure shows the 
ELGs inhabit haloes of mass between lO'^^ and IO'^A/q- The comparison with the HAM halo mass values (Table [2| shows these selections ai'e tracing the 
most massive halos at their mean redshift. 
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Figure 14. Large scale galaxy bias of the different BAO tracers ordered by 
mean redshift. The bright samples are more biased than the faint samples. 
The HOD bias is not in agreement for the faint samples and tends to be 
greater than the weak lensing estimation. 



ing discussion is therefore restricted to these tracers unless other- 
wise stated (we set aside the ELG faint samples). 



5.1 Are these selections suitable for future BAO studies ? 

The signal to noise in BAO goes as the number density (n) times 
the power spectrum, which amplitude scales as oc h- for a constant 
CTg. Thus SNRbao nP oc hb-. Thus the higher the bias (or the am- 
plitude, it is equivalent), the smaller is the number density required 
to a reach a given signal to noise, and thus the faster the survey. 
For the same signal to noise ratio, the density required using trac- 
ers with 6 = 1.5 (resp. = 2) is 2.25 (resp. 4) times smaller than 
using tracers with a bias b = \. A BAO survey using these tracers 
has thus a real advantage compared to selecting galaxies with lower 
absolute luminosity that have a lower clustering amplitude. These 
numbers constitute a zeroth order value of the galaxy bias that can 
be used in simulation to predict the accuracy of the detection and 
anticipate the requirements on geometry and density for surveys as 
SDSS-IV/eBOSS, DESspec, BigBOSS, PFS-SuMiRe, and Euclid. 
A check could be performed with n-Body simulations by trying to 
reproduce the angular clustering measurements from this paper. 

The primary role of the target selection, to select highly-biased 
galaxies, is here fulfilled. There is, however, a trade-off. The galaxy 
bias is the source of the dominating systematic on the BAO mea- 
surement; see |Mehta et al.|p01 1| > and references therein. With cur- 
rent reconstruction algorithms, and a knowledge of the bias better 
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than 10%, simulations show it is possible to reduce the system- 
atic errors down to 0.15%. This study enables to start calibrating 
simulations on observed data to understand how to cope with this 
systematic elFect. 



5.2 ELG ugri selection, an intermediate selection 

The ugri ELG color selection is an hybrid between the ugr and the 
gri selections. The clustering analysis shows the HOD of the ugri 
and ugr are similar and the value of the bias for the ugri and the gri 
are almost identical. We are thus self consistent. We checked this 
by computing the angular cross-correlation between these samples, 
and it is small on small scales, and becomes strongly correlated on 
larger scales. 



5.3 Independence of the lower redshift selections 

The redshift distribution of the CMASS, LRG-WISE and gri ELG 
demonstrates that there is a small redshift range in which the three 
samples are present 0.6 < z < 0.7; see Fig. [T] The cross-match 
between the three samples shows there is 20% overlap between the 
LRG-WISE bright sample and the CMASS and no overlap between 
the LRG-WISE bright or the CMASS and the ELG gri. The three 
samples are therefore complementary to map the large scale struc- 
tures in this redshift zone. 



5.4 Hints on luminosity, stellar mass, halo mass, HOD 
modeling and N-body simulations 

|Wake et aTjjZOl 1^ showed that stellar mass and halo mass are corre- 
lated for 1 < z < 2, and |Li et ar] ( [2006| reported that on large scales 
the galaxy clustering amplitude dependence on stellar mass is pos- 
itively correlated to the luminosity. Our selections target the most 
luminous galaxies (in their redshift range). We therefore expect to 
track high stellar masses and high halo masses. The combination 
of the clustering analysis from this work and of the estimates of 
stellar masses from Comparat et al.| j2013( ) demonstrate that these 
galaxies trace haloes with large masses (M;,^^ ~ 10"Mq) i.e. we 
comply with the general picture. 

This analysis is complementary to similar analyses based on 
faint and deep 'pencil beam' surveys, that investigate the cluster- 
ing properties of average luminosity galaxies, see |Leauthaud et aT] 
( |2012^ ; [Mostek et al.| ( |2012) . In fact our target selections aim for 
galaxies that are not in large numbers in such surveys. Given the 
trends in these papers, our results are in good agreement with what 
galaxy bias is expected for the brightest galaxies at redshift z = 1 . 

These clustering measurements are the first of its kind and are 
of prime interest to improve our knowledge of the halo occupation 
distribution of such galaxy populations. In this work, we assumed 
the angular clustering measured was representative of the cluster- 
ing of a volume limited sample at the mean redshift of the sample 
in order to derive its large scale bias. This assumption is not sim- 
ple to test with current data, in particular with current photometric 
redshifts. Though this might be one of the reasons why the HOD 
model fails to describe in details the observed angular clustering at 
small scales, in particular the satellite fractions. In fact the proper- 
ties of the selected galaxies evolve with redshift i.e. the lower red- 
shift ones are less luminous than the higher redshift ones. The HOD 
model used here does not have parameters that evolve with redshift. 



it is a fixed halo population convolved with the redshift distribu- 
tion. These considerations are oriented towards the construction of 
a mock catalog. 



6 CONCLUSION 

We measured, for the first time, the angular clustering and the cos- 
mic shear of the future samples selected for BAO studies. We have 
validated the use of this method on the CMASS sample and con- 
firmed previous measures of its galaxy bias. Moreover, we show 
that the CMASS sample is well correlated to the dark matter field. 

We demonstrate that the bright ELG samples are highly biased 
and therefore suited to detect BAO. The galaxies to be targeted for 
the near future SDSS-IV/eBOSS BAO experiment will provide a 
high signal-to-noise ratio in the galaxy power spectrum. 

This work provides a strong basis to develop simulations and 
mock catalogs of these galaxy populations in order to understand 
the details of the relation ELG - dark matter haloes. 

When pushing the method further in depth (ELG faint sam- 
ples), we are limited by the knowledge of the redshift distribution 
of the tracers and by the depth of galaxy shape catalogs for weak- 
lensing. Both methods disagree when the redshift distribution is 
poorly known, i.e., with an uncertainty of about 40%. Within the 
next few years, deeper observations and improved photometric red- 
shifts should provide the necessary data to describe with precision 
the relation between galaxies and matter probed by the faint target 
selection schemes. 
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